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We present the results of an experimental study of patterned segregation in a hori- 
zontally shaken shallow layer of a binary mixture of dry particles. As the compacity, 
C, of the mixture was increased, the evolution of three distinct phases was observed. 
We classify them as binary gas, segregation liquid and segregation crystal phases using 
macroscopic and microscopic measures. The binary gas to segregation liquid transition 
is consistent with a continuous phase transition and includes the characteristic feature 
of critical slowing down. At high compacities we observed an intriguing slow oscillatory 
state. 

1. Introduction 

In segregation, excitation via flow or shaking can, counter-intuitively, cause an initially 
homogeneous mixture of grains to de-mix [IHEI- Intriguingly, it does not always happen 
and the conditions for its occurrence are difficult to predict. Hence, a better understand- 
ing would have a major economic impact in the pharmaceutical, chemical processing and 
civil engineering industries [EJ- However, despite decades of research, a predictive model 
of the process has yet to emerge. More recently, segregation has received considerable 
attention from the physics community as an example of a challenging far from equilib- 
rium system with pattern formation Although many small scale laboratory studies 
involving vibration [ , avalanching in partially filled horizontal rotating drums [ 5^ and 
stratification in vertically poured mixtures [ IH] have been investigated and a variety of 
mechanisms proposed [ |7j , an understanding of the fundamental principles involved re- 
mains incomplete. 

In our experiments [Eld El CHI we have focused on quasi- 2- dimensional layers of binary 
mixtures of particles which are horizontally driven by the stick-slip frictional interaction 
with the surface of a horizontal tray. Vibration of an initially homogeneous mixed binary 
layer can give rise to robust segregation of the constituents, depending on the total filling 
fraction of the layer. We call this parameter the compacity. As the compacity was 
increased the evolution of distinct phases was observed. We identified them as a binary 
gas, a segregation liquid and a segregation crystal using both microscopic and macroscopic 
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measures to identify their properties. This horizontal and quasi-2-dimensional set up 
gives the practical advantage that any collective behavior is readily visualized, vertical 
compaction effects are effectively eliminated and the material is in contact with the forcing 
throughout the drive cycle. 

For this class of quasi-2D binary granular systems a qualitative segregation mechanism 
has been suggested [ [TUl HTJ [T21 ESI using the idea of excluded volume depletion as in 
colloidal systems and binary alloys [E]- This is in the spirit of the physics of complex 
systems where attempts have been made to unify descriptions of granular materials, col- 
loids, gels and foams [[T7l[TH]. We discuss the relevance of these ideas to our observations 
in some concluding remarks. 

2. The experiment 




Figure 1. (a) Schematic diagram of the experimental apparatus, (b) Trajectories of two individual 
particles over Iniin: black line corresponds to poppy seed and grey line to phosphor-bronze sphere, 
(c) Typical segregation pattern for mixture A: poppy seeds (grey regions) and phosphor-bronze spheres 
(white regions), (d) Typical segregation pattern for mixture B: sugar particles (shaded grey regions) and 
polystyrene spheres (white regions). Both patterns self-organised 15mm after starting from a homoge- 
neous mixture. 

A schematic diagram of the top view of the apparatus is presented in Fig. [D^a). It 
consisted of an accurately levelled and horizontal smooth rectangular tray, of dimensions 
(x, y) = 180 X 90mm with a flatness of less than ±5fim, on which particles were vibrated 
longitudinally. The tray was made from aluminum tool plate for increased rigidity and was 
mounted on a horizontal platform which was connected to an electro-mechanical shaker. 
Its motion was a close approximation to one- dimensional as it was constrained by four 
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Material 


Avg. Diameter {mm) 


Density {gem ^) 


Shape 


poppy seed 


1.07 (±17%) 


0.2 


flat 'kidney' shaped 


phosphor-bronze 


1.50 


8.8 


precision spherical 


polystyrene 


0.5 


1.1 


precision spherical 


sugar 


1.71 


1.6 


roughly spherical 



Table 1 

Physical properties of the particles used. 



lateral high precision linear bush bearings. The dynamic displacement and acceleration of 
the shaking bed was monitored by a Linear Displacement Variable Transformer (LVDT) 
and a piezoelectric accelerometer. The tray's acceleration was checked to be sinusoidal, 
to 0.1% in power, using the power spectrum of the accelerometer trace. 

We have investigated segregation using a variety of mixtures but here focus on results 
for two specific types: mixture A (polystyrene spheres + sugar particles) and mixture B 
(phosphor-bronze spheres + poppy seeds). In particular, we have used mixture B for the 
detailed quantitative results. The details for the particles' material properties is presented 
in Table 1. 

At this point we define the layer eompacity to be the total filling fraction of the system, 

c^Mi±3^, (1) 

xy 

where A^i and N2 are the numbers of species a and b in the layer, Ai and A2 are the two 
dimensional projected areas of the respective individual particles {Ai = (0.87±0.15)mm^ 
for poppy seeds and A2 = (1.77 ± 0.06)mm^ for phosphor-bronze spheres) and x and y 
are the longitudinal and transverse dimensions of the rectangular tray. We have chosen to 
fix the forcing at a single frequency / = 12Hz and amplitude A = ± 1.74mm and explore 
the effect of the layer eompacity. 

Throughout this study, experiments were performed in a quasi-monolayer regime so 
that the smaller particles were rarely on top of the larger ones. This was particularly the 
case for the mixtures of poppy seeds and phosphor bronze spheres since the size ratio for 
these two types of particle is ~ 2 : 3. However, some overlapping of particles did occur so 
that compacities up to C ~ 1.1 were achieved in the high packing limit. 

3. Single particle motion 

As discussed above, the forcing was strictly sinusoidal but both types of particles were 
driven by a stick-and-slip interaction with the oscillatory surface of the container. These 
interactions induced stochastic motion and trajectories of single particles with the char- 
acteristics of a 2D random walk were obtained. Typical one minute long trajectories for 
a poppy seed and a phosphor-bronze sphere are shown in Fig. ^b). The statistics of 
the motion were anisotropic with a preference for motion in the direction of the forcing. 
Note that the two types of particles responded to the driving differently, due to their 
mass, shape and surface properties. When many particles were present, collective motion 
dominated with diffusion biased in the direction of the forcing. 



4 




I , I I I I I , I 

200 tA400 600 tg 800 1000 1200 
Time (sec) 



Figure 2. (a) Time series of the average domain width, (f>, over 20min, for C = 1.057. (b)&(c) Binary 
frames of the segregation patterns (regions of poppy seeds in black and regions of phosphor-bronze spheres 
in white) at (b) tA = SOOsec and (c) ts = 702sec. 



4. Granular segregation 

In Fig. ^c) and (d) we present two examples of segregation patterns for the two mix- 
tures considered. The snapshots correspond to the segregated states which formed after 
the mixture was vibrated for a period of 15 minutes. The experimental runs were system- 
atically started from a homogeneous mixed layer using the following method. Firstly, N2 
small particles were vibrated at large amplitudes, A ~ ±5mm, creating an homogeneous 
and isotropic sub- monolayer. The large particles were then suspended above the layer, on 
a horizontal perforated plate with (57 x 28) 2mm diameter holes arranged in a triangular 
lattice. A shutter was then opened and the 1596 particles fell onto the layer of small 
particles, creating a near homogeneous binary mixture. Domains such as those shown 
in Fig. ^c) and (d) then formed. The differences in the physical properties of the two 
mixtures illustrates that, in our system, patterned segregation is a robust phenomena. 

In order to quantify the dynamics of growth of the segregated domains we focus on 
mixture B and establish two macroscopic measures: the average longitudinal width and 
the average number of domains of phosphor-bronze spheres. At this stage, the images of 
the whole layer (300 x 600 pixels) were processed such that regions of poppy seeds were 
given a value of (black) and regions of spheres a value of 1 (white). For each acquired 
frame this was calculated by scanning through each of the 300 horizontal lines and the 
length, Li, of each of the white steps was measured. We therefore define, = (L), to be 
the average value of the domain width and the average number of domains, 77, to be the 
number of such steps, per pixel line. 

A typical time evolution of </> for C = 1.057, over a period of 20min, is presented in Fig. 
|2fa). It is characterised by a fast initial segregation growth, over timescales of the order 
of ~ Imin. During this period the mechanism of the segregation process is as follows. 
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Figure 3. (a) Time-series of the average domain width, 0. (b) Time-series for number of domains, rj. 
(□) C = 0.873, (■) C = 0.729 , (•) C = 0.708 , (o) C = 0.602. Lines are fits to Eqn. ©■ AU runs were 
started from a homogeneously mixed layer. 



Initially, single large particles diffuse in a sea of the smaller ones. When two large particles 
come close together, the smaller particles cannot fit between them, and hence the pair is 
subjected to an asymmetric pressure that will keep it together. Subsequently, pairs may 
encounter others so that progressively larger clusters form. The unidirectionality of the 
driving induces an asymmetry in the segregated domains such that elongated domains of 
the larger particles (phosphor-bronze spheres or sugar particles) develop. The segregation 
domains were elongated in a direction perpendicular to the drive. The growth of the 
domains then exhibits an intermediate saturation level as can be seen in the plateau 
labelled P in Fig. Efa). 

Further coarsening of the domains can take place, albeit at longer timescales. Two 
examples of the patterns at later times are presented in Fig. Efb) and (c) at tA = 300s 
and ts = 702s, respectively. During this period coarsening has occurred via the break-up 
of one of the stripes and consequent merging with its neighbours. In contrast to the 
fast initial segregation, this slow coarsening may be associated with large scale collective 
motion in the layer of poppy seeds. 

Since there is a distinct separation of the timescales we will now focus on the fast 
initial segregation growth regime. In Fig. Efa) and (b), we present typical time-series for 
and 7], respectively, for four different values of C, over the first Amin of vibration. The 
superimposed hues are fits of, 

0(C,t) = S'^-6'^.exp (^-^^ and r/(C, t) = + 6". exp (^-^^ , (2) 

to the experimental time series of both and rj, where E"^ and are the respective 
saturation levels and and t"^, the segregation times, are timescales associated with the 
saturation of the domain growth. We denote S*^ by the segregation level. A result which 
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Figure 4. Segregation phase transition, (a) Compacity dependence of saturation levels from average 
domain width (•) and average number of domains (□). (b) Segregation time, t^. from average domain 
width (□) and average number of domains (•). (c) Inset: Segregation time from average number of 
domains zoomed around C = 0.893. Lines are guides to the eye. 



may be inferred from both sets of time sequences in Fig. El is that, when the domains 
form, there is a progressive coalescence of clusters with an increase of the average width 
and corresponding decrease in the number of domains in the system. 

5. Segregation Phase transition 

We next investigated the dependence of the saturation levels, E, and segregation times, 
tg, on the compacity of the layer for both (p and 1]. This was done by incrementally 
increasing the number of poppy seeds, A^^i, in the layer, in measured steps, while the 
number of phosphor-bronze spheres, N2 = 1596, held constant, i.e. C = C{Ni). We 
chose to perform the experiments by changing the numbers of poppy seeds. This can be 
regarded analogous to the experimental protocol in colloid-polymer mixtures of increasing 
the concentration of the polymer in solution (see Section 8 and Ref. [E])- The qualitative 
behaviour presented here was found to be robust over a range of A''2. 

The C-dependence of the segregation level is presented in Fig. EJ^a). For low com- 
pacities, up to a value, Cc ~ 0.647, the segregation level remains low and approximately 
constant around (3.25 ± 0.24)mm. This corresponds to approximately two sphere di- 
ameters and we do not classify these as domains; the system is deemed to be in a mixed 
or unsegregated state. This is essentially a collisional regime since there is enough free 
volume in the system so that particles collide and randomly diffuse, without aggregation. 
We denote this regime by a binary gas. 

As C is increased, past Cc, clusters of the larger particles form mobile segregation 
domains, of increasing 0, emerge in a nonlinear way. Above Cc, the clusters of larger 
particles were mobile and their movement was reminiscent of oil drops on water as they 
flowed, merged and split and the motion of the particles within the clusters was highly 
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agitated. As C was incrementally increased there was a gradual decrease of the overall 
motion of segregation domains so that they become less mobile and increasingly more 
stable. Also, for larger values of C, there was a significant decrease of the agitation of the 
spheres within the domains. 

The average number of domains, shown in Fig. |3 remains approximately constant up 
to the same value Cc, after which it rapidly decreases as C is increased. This is expected 
for increasingly wider domains since the number of phosphor-bronze spheres is constant. 
This scenario is consistent with the existence of a continuous phase transition for granular 
segregation of the binary mixture, as the layer compacity is incrementally increased, with 
a critical point for segregation, Cc — 0.647. 

In Fig. El^b) we present the segregation times obtained from the fits of Eq. Q to 
and ri{t) as a function of C. Note that we only define for C > Cc, since this 
is the parameter range over which segregation occurs. As expected for a continuous 
phase transition, tg diverges, as C is decreased from above, near Cc, a phenomena usually 
referred to as critical slowing do'wn['^^. This adds further evidence for the existence of a 
phase transition. 

As the compacity is increased, past Cc, this slowing down of the segregation persists 
up to C = 0.893. At this point there is a minimum in the segregation time, as shown in 
Fig. O^c), with tg = 11.41s. For greater values of C there is a qualitative change in the 
segregation behaviour and the layer takes increasingly longer to segregate. At this large 
values of C, stable crystalline stripes form which are separated by dense regions of poppy 
seeds. 

The above scenario suggests that the segregated phase, for C > Cc, can be sub-divided 
into two further regimes: one in which the speed of segregation increases with C (0.647 < 
C < 0.893) and another in which the opposite occurs (C > 0.893). We denote the first 
by segregation liquid and the second by segregation crystal. 

6. Binary gas, segregation liquid and segregation crystal phases 

To further justify this idea of three segregation phases we turn to the structural infor- 
mation given by the radial pair distribution function for the large spheres, g{r). 



where rij is the separation of the i and jth particles and the angled brackets denote a time 
average. A{r) is a normalization constant such that g{r) = 1 for a uniform distribution 
of particles. The data for the calculation of g{r) was obtained 3min after the driving was 
switched on. Time averaging was then performed over Imin, for 1500 temporal realisations 
of the layer's configuration. 

In Fig. Efa) we show three curves of g{r) where each represents the segregation phases 
discussed above. Typical configurations of the spheres, for each of the phases, are pre- 
sented in the zoomed photographs in Fig. El^b),(c) and (d). Note that the photographs 
have been image processed such that only the phosphor-bronze spheres are visible. For 
the binary gas (C < 0.647), g{r) has a peak a.t r/d = 1 which quickly decays at large 
distances, as expected for a disordered gas. Increasing C results in a monotonic increase of 




(3) 



i j^i 
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Figure 5. (a) Radial distribution functions, g{r), for the phosphor-bronze spheres: (a) binary gas 
(C=0.503), (■) segregation hquid (C=0.708) and (•) segregation crystal (C=1.037). The vertical dashed 
lines correspond to the expected location of the peaks for an hexagonal lattice in 2D. (b) Zoomed binary 
gas. (c) Zoomed segregation liquid, (d) Zoomed segregation crystal. The photographs have been image 
processed such that only the phosphor-bronze spheres are visible. 



the height of the first peak, g{d). At intermediate compacities (0.647 < C < 0.893), liquid- 
like behavior is observed, with g{r) peaked at 1, 2 and 3 particle diameters. The positions 
of neighbouring spheres are correlated and the maxima may be associated with shells of 
neighbours. The oscillations are rapidly damped, showing the decay of short-range order. 
This behaviour is commonly seen in hard sphere liquids and was first observed experi- 
mentally by Bernal [CSl- In the segregation crystal phase (C > 0.893) two further peaks 
emerge near r/d = \/3 and r/d = a/7, characteristic of a hexagonally packed crystal, in 
two dimensions. At these high compacities, individual spheres in the segregation domains 
rattle inside the cage formed by its neighbours. 

7. Oscillatory states 

At the highest values of the compacity where robust crystalline stripes form, we have 
observed intriguing well defined oscillations of the segregation patterns. The phosphor- 
bronze stripes bend, periodically, backwards and forwards in a reproducible way. For 
each of the regions of poppy seeds in between two phosphor-bronze stripes, there is a 
complementary collective motion in the form of a single large vortex. The direction of 
rotation of the vortex is reversed every half-cycle of the oscillation. In Fig. inta)-(c) we 
show three photographs of the granular layer which were captured 40 min after having 
started from a homogenous mixture and correspond to the beginning of the cycle, half- 
cycle and full cycle points of one oscillation, respectively. 

We illustrate the periodic nature of the oscillations, in Fig. Efb) with a space time 
diagram of the process. This was constructed by sampling along a single line in the x- 
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Time At={1.5 hrs) 



Figure 6. Photographs of the granular layer showing the oscillatory bending of the 
phosphor-bronze stripes (regions in white), at different times: (a) beginning of cycle, (b) 
half-cycle, (c) end of cycle, (d) Space time diagram constructed from stacks of the cut 
line, L in (c), over 1.5/irs. 



dimension, and progressively stacking the lines over a period of 1.5hrs. The sampling 
line was positioned at 1/3 of the y-length of the tray as indicated by the solid white line 
superposed in Fig. Efc). A well defined cycle is clearly observed and it has a period 
of T = 286s. This suggests the existence of a strong coupling between the oscillation 
of the stripes through bending and the vortex inversion in the poppy seed's layer. It is 
interesting to note that the timescale associated with these oscillations is ~ 3500 times 
slower than the driving (at 12Hz) excluding the possibility of a simple resonance in the 
driving/granular-layer system. Hence, the mechanism which gives rise to this oscillatory 
phenomena remains, as yet, unexplained. 

8. Discussion and conclusion 

It is interesting to note that there are some analogies between segregation in our binary 
granular system and aggregation in binary colloids [132] and emulsions ['^. In particular, 
when the size ratio between the large colloidal particles and the radius of gyration of the 
polymer is larger than ~ 0.3, colloid-polymer mixtures display colloidal gas, liquid and 
crystal phases of the large spheres [ |2I1 122] • Aggregation in such systems is usually ex- 
plained using an excluded volume depletion argument introduced by Asakura and Oosawa 
[ I2S1 in the context of binary hard spheres. This entropic argument is equivalent to the 
mechanistic view that if two large particles are close enough so that no other particle (or 
polymer) may fit between them, they will be subjected to an asymmetric osmotic pressure 
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that leads to an effective attractive inter-particle force. In colloid-polymer mixtures the 
strength of the interaction can be tuned by changing the concentration of the polymer in 
solution, therefore inducing gas-to-liquid and liquid-to-crystal transitions. These resemble 
the granular phases we have observed by incremental increase of the compacity. However, 
an important difference is that our granular system is both driven and dissipativc and 
therefore far from equilibrium. In addition, Brownian motion is an intrinsic part of the 
equilibrium dynamics of colloidal systems. Hence a direct connection between the two 
systems remains an open question. 

In conclusion, our results provide quantitative evidence for the existence of three phases 
- binary gas, segregation liquid and segregation crystal - in granular segregation of hori- 
zontally excited binary mixtures. At high compacitics, intriguing regular oscillatory states 
were observed. We show that segregation undergoes a continuous phase transition that 
occurs for compacities above a critical value and exhibits critical slowing down. We have 
presented both macroscopic and microscopic measures which are self-consistent. Analo- 
gies of this phase behaviour may be drawn with other depletion driven self- assembling 
binary systems. This raises the possibility of importing ideas from binary colloids in equi- 
librium to formulate new models for granular segregation. 
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